> TGRS-2007-00582 < 


1 


Differential Radiometers Using Fabry-Perot 
Interferometric Technique for Remote Sensing 

of Greenhouse Gases 


Elena M. Georgieva, William S. Heaps, Emily L. Wilson 


Abstract — A new type of remote sensing radiometer based 
upon the Fabry-Perot interferometric technique has been 
developed at NASA’s Goddard Space Flight Center and tested 
from both ground and aircraft platforms. The sensor uses 
direct or reflected sunlight and has channels for measuring 
column concentration of carbon dioxide at 1570 nm, oxygen 
lines sensitive to pressure and temperature at 762 and 768 nm, 
and water vapor (940 nm). A solid Fabry-Perot etalon is used 
as a tunable narrow bandpass filter to restrict the 
measurement to the gas of interest’s absorption bands. By 
adjusting the temperature of the etalon, which changes the 
index of refraction of its material, the transmission fringes can 
be brought into nearly exact correspondence with absorption 
lines of the particular species. With this alignment between 
absorption lines and fringes, changes in the amount of a species 
in the atmosphere strongly affect the amount of light 
transmitted by the etalon and can be related to gas 
concentration. The technique is applicable to different 
chemical species. We have performed simulations and 
instrument design studies for CH 4 , 13 C0 2 isotope, and CO 
detection. 

Index Terms — Absorbing media, Atmospheric measurements, 
Fabry-Perot interferometers, Optical interferometry, Remote 
sensing. 


I. Introduction 

P recise detection of carbon dioxide in the atmosphere is 
of great interest due to its impact on trapping the long 
wavelength radiation emitted from the Earth’s surface. 
There are different scenarios about the detailed causes for 
the global warming of the Earth’s atmosphere. The net 
impact of anthropogenic aerosol emissions and greenhouse 
gases is a complex study [1]. Recent climate models focus 
more on the effects of aerosol, which can absorb solar 
energy and warm the atmosphere. However the indirect 
effects also occur when aerosol change the cloud optical 
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properties and reflect the solar energy back into space 
[2, 3, 4, 5]. Studies have shown that aerosol and greenhouse 
gas forcing are almost of the same order of magnitude [1,3]. 
However while the anthropogenic aerosols have short 
lifetime, the carbon dioxide is a long-lived constituent. The 
concentrations of the principal anthropogenic greenhouse 
gases (C0 2 , CH 4 , N 2 0) have increased substantially during 
the industrial period. The carbon dioxide concentration has 
increased by more than 95 ppmv in the last 150 years due to 
the emissions of the modem industrial revolution and the 
mixing ratio has become 380 ppmv [6, 7]. In this same time 
frame, the concentration of atmospheric methane has 
increased by 150 percent from approximately 700 to 1,745 
ppbv [8, 9]. Although found in lower concentrations than 
carbon dioxide, methane is about 20 times more effective in 
absorbing infrared radiation. While the majority of C0 2 
variability occurs in the lower atmosphere (-1000 to 800 
mbar), satellite instruments measure the total atmospheric 
column. Since sources and sinks at the surface represent a 
small perturbation to the total column, a precision of better 
than 1% is required [10, 11]. That number comes from the 
transport inversion model experiments, which indicate that 
global column measurements with a precision better than 
1% (3ppmv on the 380 ppmv background) on a time scale 
of one month is the science requirement to improve surface 
flux estimates beyond the capability of the existing network 
[11]. The natural geographic distribution and temporal 
variability of C0 2 sources and sinks are still not well 
understood, and more monitoring instruments are needed to 
quantify the carbon dioxide dynamics and to help predict 
climate change [12,13]. The water vapor concentration 
resulting from evaporation of water from land and ocean is 
enhanced by carbon dioxide. As carbon dioxide absorbs 
radiative heat from the Earth’s surface - warming the 
atmosphere, atmospheric water vapor increases because the 
warmer air can hold more water [14]. The increased water 
vapor leads to an increase in the greenhouse effect and 
further increase in temperature; the cycle is repeated until 
equilibrium. For climate research, high precision 

(error<5%) water vapor measurements are needed with 
global coverage [15]. 

Our knowledge of atmospheric carbon processes comes 
from satellites (Atmospheric Infrared Sounder, 
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SCIAMACHY onboard ENVISAT satellite [16]) and 
ground network measurements. The latter consists of two 
main sources- the long term, in situ C0 2 measurement 
program led by the NOAA Climate Monitoring and 
Diagnostics Laboratory and the TransCom 3 transport/flux 
estimation experiment [17, 18]. These in situ measurements 
are very precise (uncertainties on the order of 0. 1 ppm) and 
accurate, but are necessarily limited in time and space. 
Rayner and O’Brien [11] predicted that the required 
accuracy for monthly averaged C0 2 column data needs to 
be better than 2.5 ppmv on 8° x 10° footprint for 
comparable performance to a moderate surface network. 
Such a precision is difficult to achieve. 

The instrument that we have developed and tested on two 
flight campaigns detects the absorption of C0 2 , 0 2 and H 2 0 
gases using direct or reflected sunlight [19, 20, 21, 22, 23]. 
The sensor makes use of two features of the Fabry-Perot 
interferometer to achieve high sensitivity: high spectral 
resolution resulting from matching the width of an 
atmospheric absorption feature to the instrumental band 
pass, and high optical throughput enabled by simultaneous 
use of multiple spectral lines. For any species that one 
wishes to measure, this first feature is available while the 
use of multiple spectral features can be employed only for 
species with suitable spectra and freedom from interfering 
species in the same wavelength region. Throughput of an 
interferometer can be as much as two orders of magnitude 
larger than a spectrometer with a grating with the same 
resolution, which means larger signals or better signal-to- 
noise ratio (S/N). Our sensor records data every 0.1 s, has 
high spectral selectivity and has only one moving part 
(chopper), which makes it an ideal candidate for an airborne 
or satellite mission. It can work also as a ground based 
sensor using direct light from the Sun. The current goal is to 
develop a precise, inexpensive, ground based device 
sensitive enough to measure column averaged C0 2 and 
suitable for wide deployment as a validation instrument for 
the Orbiting Carbon Observatory (OCO) satellite [24, 25, 
26, 27]. This will fill the gaps of the sparse ground network 
and help to understand the global distribution of C0 2 and its 
role in climate change. Parallel to that we are developing a 
new sensor for methane detection and working on 
simulations for systems to measure CO and 13 C0 2 isotope. 


aperture. The incident light in the C0 2 channel is pre- 
filtered at a central wavelength of 1570 nm. The light then 
splits between the Fabry-Perot and Reference sub-channels 
with 90% of the light going to the Fabry-Perot sub-channel. 
In the Fabry-Perot sub-channel the light passes through the 
Fabry-Perot etalon mounted in a temperature controlled 
oven for fine FSR (free spectral range) tuning. By adjusting 
the temperature of the etalon, which changes the index of 
refraction of its material, the transmission fringes can be 
brought into nearly exact correspondence with absorption 
lines of the particular species (Fig.2). A set of off axis 
parabolic (OAP) gold mirrors focuses light onto InGaAs 
detectors for the C0 2 sub-channels and Si detectors for the 
two oxygen channels. For bringing light into the instrument 
we use fiber bundles, which intentionally scramble the 
incident light to remove spatial information. The Fabry- 
Perot etalons in this instrument are solid fused silica etalons 
and they have different thickness for each instrument 
channel. For the C0 2 channel the etalon has a free spectral 
range (FSR) of 0.306 nm, a refractive index of 1.443 at 
L=1571 nm, and a clear aperture of 50 mm with highly 
reflective coatings on both surfaces. The light that passes 
through the etalon undergoes multiple reflections on each 
inside surface, creating an interference pattern of equidistant 
fringes as a function of frequency. The width of the pass 
bands depends on the coatings quality and on the flatness 
and parallelism of the surfaces. The Fresnel formalism gives 
the reflected and transmitted amplitude components for the 
light wave as functions of the optical constants of the two 
media and the angle of incidence. The resulting intensity 
follows the Airy function distribution pattern. The ideal 
Fabry-Perot etalon with perfectly flat surfaces transmits a 
narrow spectral band and the energy transmission 
coefficient I T is given by [28, 29]: 


l T 


(l-R) 2 


i + - 


4 R 

(1 - R ) 2 


2 ( iTtndcosO 


( 1 ) 


where X is the wavelength, n is the refractive index, d is the 
thickness of the etalon and 0 is the angle of incidence within 
the cavity, T is the intensity transmission coefficient for 
each coating and R is the intensity reflection coefficient. 
The transmission is periodic. 


II. DESCRIPTION OF INSTRUMENT AND BASIC 
THEORY 

The incoming light is split into the three wavelength 
bands with a dichroic beam splitter and directed into C0 2 , 
0 2 pressure-sensitive and 0 2 temperature-sensitive 
channels. The optical schematic for the C0 2 channel of the 
FP instrument is shown in Fig.l. In this design, incoming 
light is focused onto a 2 mm diameter aperture to define the 
field of view. The light is modulated at -400 Hz with a 
chopper and then re-collimated as it emerges from the 


III. EXPERIMENTAL 

A. Data taken with the flight version of the FP 
radiometer 

Carbon dioxide and oxygen column measurements using 
FP radiometer have been demonstrated in laboratory, 
ground based and airborne experiments. Water vapor 
channel was added recently and has been tested in the lab 
and ground based measurements. The FP radiometer is 
taking daily measurements at Goddard as a ground based 
sensor and the C0 2 , 0 2 and H 2 0 columns are measured 
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through absorption of light by those species in the 
atmosphere directly between the sun and the instrument. We 
accomplished this by collecting light with a small telescope 
fixed to an equatorial mount, aligned to track the sun 
throughout the day. An optical fiber coupled at the rear of 
the collimator brings light into the instrument. The initial 
tests indicate that the instrument can detect changes in the 
C0 2 column as small as 2.3 ppm with a one second average 
and better than 1 ppm in less than 10 seconds of averaging. 
The precision of the oxygen column pressure changes is as 
small as 0.88 mbar. Water vapor precision is 2 percent for 
one second of averaging. Figure 3 illustrates a typical data 
set collected over the course of a clear day with the 
instrument. The ratio of Fabry-Perot to reference signals is 
inversely proportional to absorption for C0 2 and H 2 0 and 
direcdy proportional for 0 2 which is plotted on an inverted 
scale. Figure 4 shows the same data plotted as a function of 
the calculated air mass (the pathlength of the sunlight 
through the atmosphere between the sun and the 
instrument). The air mass decreases all morning until local 
noon and then increases until sunset. The air mass m is 
expressed by the Chapman function in the IDL program 
(Interactive Data Language) for data analysis so we can 
account for the sphericity of the Earth. 



where m is the air mass, x is the optical depth, a is the 
extinction coefficient as a function of altitude z, 0 is the 
zenith angle of the sun, n is the atmospheric index of 
refraction. For very small solar zenith angles the air mass is 
equal to the secant of the zenith angle and that is the plane- 
parallel approximation. 

The 0 2 data essentially lie on top of one another because 
the atmospheric pressure (and hence the 0 2 column) was 
constant on this day. The C0 2 channel shown in blue has 
some slight variability. C0 2 is expected to change very 
slightly on a diurnal basis because of the consumption of 
C0 2 by plants during the day. C0 2 has also a seasonal 
variability with maxima in spring and minima in the late 
summer [30] .Those variations reflect changes in the 
biological sinks and can range from 1 ppmv in the Southern 
hemisphere (SH) to more than 20 ppmv in the Northern 
hemisphere (NH). 

The variation in the water vapor column is shown to be 
real in a comparison with total column measurements by 
NOAA using GPS (Global Positioning System) technology. 
This result verifies that instruments of this type function 
properly and can be calibrated precisely. Figure 5 shows 
this inter-comparison of H 2 0 measurement by the Fabry- 
Perot instrument system and by a NOAA GPS based sensor 
located in Annapolis, Maryland for 5 days in late June and 


early July. The FP radiometer only operates during the day. 
Agreement is consistent with variability expected for sites 
-30 miles apart. 

For the airborne or satellite measurements the light 
passing through the atmosphere reflects on the Earth’s 
surface before entering the instrument platform. The flight 
hardened version of the instrument was tested at two flight 
campaigns at NASA Dryden Research Center and at New 
Hampshire for the Polar Aura Validation campaign on 
NASA’s DC-8 research airplane. Flights were conducted 
over a variety of surfaces (vegetation, water, snow) and 
under different atmospheric conditions. An in situ 
instrument provided C0 2 profiles up to the flight altitude. 
Temperature and density profiles were available from the 
aircraft data system. Comparison of FP radiometer data at 
different altitudes with the integrated in situ profiles formed 
the basic measure of the remote sensing sensitivity and 
provided a reference calibration. Data for the airborne tests 
were collected viewing in nadir direction. Simultaneous 
measurements of carbon dioxide and oxygen are shown in 
Figure 6(a,b). In both plots, the ratio of Fabry-Perot to 
Reference signals are compared to the altitude of the plane, 
corresponding to changes in the total path length of sunlight 
through the absorber. For C0 2 channel the ratio of Fabry- 
Perot to Reference signals (blue) are inversely proportional 
to the column of C0 2 measured and consequently the 
altitude from which measurements are made (red). As 
expected, changes in ratio clearly track changes in altitude. 
We have used the change in ratio with altitude to estimate 
the instmment response. We obtained a value of -0.44 for 
the change in ratio divided by the change in total column - 
indicating a ratio change of approximately 40% for a 
change of one air mass. Improvements to the instrument 
since this time have increased this value. Results from the 
flight experiments were reported [19, 20, 22]. The airborne 
instrument using light reflected off the ground has a 2% 
precision. The precision for the airborne experiments is 
reduced because the aerosol and cirrus scattering 
significantly alter the sunlight path through the atmosphere 
and hence the spectral attenuation. Atmospheric scattering 
can produce path length changes much larger than 1% and 
depending on the distribution of scattering particles and the 
elevation of the sun the changes can be either positive or 
negative. The OCO instrument under development at Jet 
Propulsion Laboratory (JPL) will measure the full spectrum 
of C0 2 and 0 2 and the data are expected to provide enough 
information to permit a correction for the atmospheric 
scattering effects. Because we use a different technique the 
approach for our instrument is to employ the glint off water 
for the airborne measurements. Glint has two advantages - 
it is bright and therefore will dominate the weak scattering 
from aerosols and thin cirrus. Secondly the path length is 
known very precisely since the viewing angle to the glint 
can be measured and the location of the Sun can be 
calculated. Theoretical simulations show that our 
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instrument, employing the glint is capable of measuring the 
C0 2 column with a precision better than 0.3% for space or 
airborne platforms. 

B. New sensor design 

We are in a process of developing a reduced size and 
volume portable channel for C0 2 capable of doing 
simultaneous measurements with the AERONET (Aerosol 
Robotic Network) at NASA, Goddard. Fabry-Perot 
radiometer as a ground-based sensor has the potential to fill- 
in the sparse network of surface measurements - key in 
balancing the carbon budget. The instrument’s low cost will 
enable its wide implementation in order to provide a 
broader range of intercomparison sites for the OCO 
instruments. For that purpose we designed a radiometer 
(Fig.7) which has a simpler configuration than the tested 
airborne instrument (Fig.l). The small radiometer consists 
of two sub-channels - the sub-channel that includes very 
narrow bandpass filter and the sub-channel, which detects 
the total light intensity limited by the wide bandpass filter 
spectral range. Laser scans of the two filters are shown on 
Fig.8. The C0 2 absorption lines from HITRAN database are 
also shown. The narrow filter is fitting into the gap of the P 
and R branches of the C0 2 absorption band near 1.57 
micron. All other C0 2 absorption lines are blocked by the 
narrow filter this way. They are however detected by the 
wide filter making it the active channel for this sensor. 
Sensitivity of the ground-network version of the field 
instrument is expected to be comparable to that of the 
previous ground instrument. In Table 1 we summarized 
some of the instrument parameters. 

C. Performance Simulations and Expected Results for 

Methane CH 4 and carbon isotope 13 CO 2 sensors 

A sensor for methane is under development. We modeled 
the expected performance of the instrument (Fig.9, top) as 
follows: First, we generated a synthetic spectrum for 
reflected light containing spectral features for C0 2 , H 2 0, 
and CH 4 from the HITRAN database [31]. Initially we 
assumed an average column mixing ratio of 1.8 ppm for 
methane. We calculated the solar flux signal received in the 
two subchannels (the FP and the Reference) for light that 
has passed through the atmosphere and been collected by 
the suntracker/telescope. After perturbing the total methane 
column by 1% we recalculate the signals. We compare the 
change in the ratio (FP/REF) caused by the perturbation to 
the detector noise of this ratio. The free parameters of the 
instrument design (eg. FP free spectral range, FP finesse, 
filter bandpass, filter edge slopes, instrument FOV) are 
varied in an effort to minimize response to changes in C0 2 
and H 2 0 while maximizing response to CH 4 . In the 1.6 
micron region CH 4 , C0 2 , carbon monoxide (CO) and H 2 0 
have significant absorption features. Our radiometer design 
is based on the attempt of increasing the signal for CH 4 and 
reducing that for others interfering species. It will have the 
capability of measuring CH 4 in both the Earth and Mars 
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environment. For CH 4 measurements on Mars the thickness 
of the solid etalon is calculated to be 0.15 mm, the free 
spectral range (FSR) is 6.6 nm and the prefilter will be 1 1.4 
nm wide. For the Earth measurements the solid etalon will 
have 0.359 mm thickness, FSR=2.6 nm and the prefilter will 
be 1 1.5 nm wide. On Earth only water vapor has significant 
abundance in the same spectral region. For a column 
average methane of 1.8 ppm and of 5000 ppm for water 
vapor this instrument has a theoretical signal to noise ratio 
of 1200:1 for methane and less than 1:1 for water in a one 
second averaging period. This means if the water vapor 
doubled from 5000 ppm to 10000 ppm it would only cause 
the instrument to make an error less than that arising from 
detector noise. Simulations for 13 C0 2 and CO were done in 
the same way. 

IV. CONCLUSION 

We presented here a prototype passive Fabry-Perot based 
instrument for absorption measurements of total column 
carbon dioxide using sunlight at 1.58(xm. This instrument 
has been demonstrated as a ground based sensor when used 
in conjunction with a positioning system to follow sun 
movement. Signals from these data trials were found to 
consistently track changes in air mass. Additionally, 
laboratory experiments and airborne measurements validate 
the instrument’s sensitivity to changes in C02. Additional 
water vapor channel was added to the system and prove to 
give reasonable results. New channels for carbon monoxide 
(CO), methane (CH4) and carbon isotope ( 13 C0 2 ) are under 
development. The theoretical simulations for them show 
that sufficient signal-to-noise ratios can be achieved. 

The Fabry-Perot radiometer is a mature technology for 
ground based or airborne measurements capable of 
providing high precision and high spatial resolution data of 
carbon dioxide, oxygen, and water vapor. It is doing rapid 
measurements of the entire atmospheric column, has a 
relatively simple design and can use both nadir and glint 
viewing modes. The measurements can be used directly to 
validate products from a number of satellites currently 
operated by NASA and ESA. They can also be used as 
additional checks on the measurements of some of the lidar 
and in situ sensors. 
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TABLE 1 

Instrument parameters 

The optical aperture for the airborne FP radiometer (Fig. 1) is 0.051 m. For 
the small version FP (Fig.8) the size is reduced to 0.025m. The FOV is 
0.086 deg for both instruments. The dimensions of the airborne FP 
radiometer are 0.72 x 0.35 x 0.35m and the dimensions for the small FP 
radiometer are 0.27 x 0.2 x 0.09 m. 



Bandpass 

filters 

FSR of the 
FP etalons 

Resolution 

AAA 

Sensitivity 

Channels 

FLIGHT VERSION of the FP RADIOMETER 

Carbon 

dioxide 

(C0 2 ) 

1567- 

1574nm 

0.306 nm 

13898 

2.3ppm in 1 
sec (direct 
light), 2% 
reflected light 

Oxygen 

(0 2 ) 

pressure 

760-764 

nm 

2.212 nm 

2422 

0.88 mbar 

Oxygen 

(0 2 ) 

temperature 

767-77 lnm 

0.575 nm 

4366 

0.4 deg 


SMALL VERSION of the FP 
RADIOMETER 


Water 

vapor 

(H 2 0) 

917.7- 
llOlnm & 
925-970nm 

Bandpass 
filter, 
holograph 
ic notch 
filter 


2% for 1 sec 

Carbon 

dioxide 

(C0 2 ) 

1568- 
1584nm 
1573.5- 
1578 nm 

Bandpass 

filters 


Under test 

Methane 

(CH 4 ) 

1634- 

1644nm 

2.6 nm 

6559 

Under test 

Carbon 

isotope 

( 13 CC>2) 

2111- 

2121nm 

0.4227 

nm 

42000 

Under test 
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Fig. 1. Optical and electronic schematic for the CO 2 channel of the Fabry- 
Perot instrument. Incoming light is collimated by two off axis parabolic 
gold mirrors (OAP’s). The light is divided into Fabry-Perot (detector 2) 
and reference (detector 1) sub-channels, measuring changes in CO 2 
absorption and solar flux, respectively. The ratio of these signals can then 
be related to changes in the CO 2 column. 



Wavelength, nm 


— F-P channel 



Wavelength, nm 


Fig. 2. Alignment of etalon transmission bands (blue) with CO 2 absorption 
lines (pink) can be adjusted by changing the etalon temperature. An etalon 
temperature of 54°C (top) shows better alignment than at 42°C (bottom). 


Intensity reference channel, V Intensity reference channel, V 
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WATER VAPOR RATIO 



6 t 10 12 14 il ih M 


LOCAL DAYTIME 

Fig. 3. The FP (Fabry-Perot) to Reference channel ratio is plotted as a 
function of local time for an O 2 instrument (yellow), a CO 2 instrument 
(blue), and a water vapor instrument (red).The principal variation is that 
arising from the change in air mass as the sun rises and sets throughout the 
day. Mesurement is done June 30, 2006, data were collected every 0. 1 sec 
and averaged for every 120seconds. 


WATER VAPOR RATIO 



AIRMASS 


Fig. 4. Here the ratios from the previous figure are plotted against the 
calculated air mass. The AM and PM values for O 2 are virtually the same 
and cannot be distinguished. The CO 2 channel shows some slight 
variation throughout the day probably arising from consumption by plants. 
The water vapor channel shows considerable variation consistent with 
typical changes in humidity throughout a day. 



Fig. 5. Water vapor data taken June 29, 2006 at Goddard using FPICC 
(blue dots) compared to Annapolis water vapor data (red dots). The FPICC 
data and the NOAA data using GPS technology show remarkable 
agreement. 


3) Results from carbon dioxide channel 
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pj Results from oxygen pressure sensing channel 
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Fig. 6. The instruments response to the altitude in meters at airborne 
mission, a) The CO 2 ratio of Fabry-Perot to Reference signals is shown 
compared to changes in altitude on a May 14, 2004 flight. As the aircraft 
descends the photon pathlength gets shorter so the absorption in the Fabry- 
Perot subchannel decreases, b) The Ratio for the pressure-sensing channel 
(O 2 ) and air mass as a function of time are shown. The embedded graph 
shows the blown images of the measured Ratio and altitude for four 
minutes. The standard deviation is calculated to be 0.00393. 


Oxygen pressure channel ratio 
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Fig. 7 . A representation of the physical design of the small-size CO2 
instrument. 
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Fig. 8. Laser scan of the two filters for the small CO2 radiometer. The 
C 02 absorption lines from HIITRAN database are also shown. The narrow 
band versus wide band sub-channels ratio will be sensitive to CO2 changes 
in the atmosphere. 


Methane Absorption Lines with Fabry-Perot Transmission Lines and Prefilter Shape Overlay 



WAVELENGTH 

Fig. 9 . The transmission bands of the solid Fabry-Perot etalons (yellow 
and green) and pass bands of the narrow prefilters (red) are shown. The 
top picture is representing the model simulations for the Earth atmosphere 
and the bottom picture - for the Mars atmosphere. The Martian 
atmosphere is less than 1 % as dense as Earths. The absorption features of 
CH4 are correspondingly weaker. Fabry-Perot transmission bands are 
shown aligned with some CH4 lines. There is also a strong absorption from 
water vapor (H2O) which can be cancelled by positioning of the etalon and 
the prefilter. 


Line Strength 





